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Abstract 

QJ , We study the scattering theory for the Gross-Neveu model on the half-line. We 



find the reflection matrices for the elementary fermions, and by fusion we compute 
the ones for the two-particle bound-states, showing that they satisfy non-trivial 
bootstrap consistency conditions. We also compute more general reflection matrices 
for the Gross-Neveu model and the nonlinear sigma model, and argue that they 
correspond to the integrable boundary conditions we identified in our previous 
paper M. 
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1 Introduction 



After the seminal paper by Ghoshal and Zamolodchikov Q (see also ||) on integrable 
quantum field theories with boundary, a lot of work has been done extending it, especially 
in the analysis of different models. In this paper we study one of these extensions, namely 
the O(N) Gross- Neveu (GN) model |§ on the half-line, which is closely related to the 
boundary O(N) non-linear sigma (nkr) model studied by Ghoshal in [|J. Recently [Q 
we have found new integrable boundary conditions (bc's) for the GN and nlcx models, 
based on the microscopic (lagrangian) description of these models. Here we find general 
(diagonal) solutions for the boundary Yang-Baxter equation and propose a one-to-one 
correspondence between these solutions and the boundary conditions we found in ||. 

This paper is organized as follows. In the next section we briefly review the GN 
and nlcr models, and their exact S'-matrices; in section 3 we find the exact reflection 
matrices for the GN model with the physically interesting free (Neumann) and fixed 
boundary conditions; in section 4 we solve the boundary Yang-Baxter equation for the 
more general boundary conditions we found in our previous paper ||, and establish the 
correspondence between those boundary conditions and the reflection matrices we find; 
in the final section we present our conclusions. We also have an appendix where we write 
the exact amplitudes for the scattering between bound-states of elementary fermions in 
the GN model [0, ||. Their computation is a simple exercise in fusion, but since we need 
their explicit form in this paper, we give the results in the appendix. 

2 The Gross-Neveu Model 

The GN model is defined by the following lagrangian 

£ 9 „ = ^# + y(^) 2 , (2.1) 

where ip is a iV-component massless Majorana fermion in the fundamental representation 
of O(N). In the above equation ipip should be understood as J2iL% f/'V* and so on. In 
this paper we will refer to the O(N) GN model simply as GN model. It is useful to write 
the GN model lagrangian in light-cone coordinates, where it reads 

C gn = 2^\id^\ + 2^1*0+^1 + 2g 2 (tfj i + ^ i _) 2 . (2.2) 

The particle spectrum of the GN model || is composed by the O(N) vector multiplet 
of elementary particles included in the lagrangian (the "elementary fermions" ) and a set 
of O(N) multiplets (scalar and higher rank antisymmetric tensors) of increasing mass, 
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which can be thought of as bound-states of a number of elementary fermionsQ. In what 
follows we will restrict our considerations to the sector of the theory containing only 
the elementary fermions (denoted by the symbols Ai, i — 1, . . . , N) and the two-fermion 
bound-states: the isoscalar particle (denoted by B), corresponding to the bound-state 
in the isoscalar channel of the S'-matrixof elementary fermions, and the antisymmet- 
ric multiplet (denoted by B^), corresponding to the bound-state in the antisymmetric 
channel. 

The exact S-matrixfor the elementary fermions was found by Zamolodchikov and 
Zamolodchikov in and we quote it here for further reference. The Faddeev- Zamolod- 
chikov algebra is 

n 

it=i 

+ (T2{e i2 )M^)M^) + ^3(012)^(02)^(00 , (2.3) 

where 9 12 = 9\ — 9 2 ■ The <Ji{9) are given by 

vM) = — i^-MQ) > °*lfi) = ~%(9) , (2.4) 

ITT — (7 (7 

and 

sinhtf + isinA 

<72{9) = —— r— - &Z{6) , 2.5 

smh 9 — i sin A 

with 

rf A _ it) pc I -it) pc I + A + iL) pel + M.) 

2l) r(i + A-f)r(-f)r(i + A + i|)r ( i + f)' 

where we introduced A = 2n/(N — 2). The pole at 9 = iX in the CDD prefactor in 
equation |275] gives rise, by bootstrap, to the full exact spectrum of bound-states. 

The S"-matrixelements for B and Bij particles can be obtained from the elementary 
fermion 5-matrixby fusion, using the identity 



Ai{9 - i\/2) Aj(9 + iA/2) = Sij B{9) + ^{N-A)B ij {6) . (2.7) 

A partial list of these amplitudes can be found in the appendix. 

The GN model is closely related to the O(N) nla model; they share very similar 
properties, and in particular the elementary particles in both models have the same S- 
matrix, up to a CDD factor. The O(N) nlcr model is defined by the following lagrangian 

Cni* = TT-d^ii ■ d»n , (2.8) 

where n is a vector in ^-dimensional space, subject to the constraint n • n = 1. The 

exact S'-matrixfor the elementary particles in the nlcr model is given by equations |2.3| - |2.6| 

with o 2 substituted by a\ flit)] , and consequently there are no bound-states in this model. 

1 In fact, the spectrum contains also kink states, associated to the spontaneous breaking of the chiral 
symmetry JoJ . 
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3 Reflection Matrices 



In this section we will compute the reflection matrices for the GN model with free and 
fixed boundary conditions, following Ghoshal's analysis for the nlcx model ||. Before we 
proceed let us note that what Ghoshal means by fixed boundary condition is not what 
we meant by Dirichlet boundary condition in our previous paper || (and in this one). 
Ghoshal's condition corresponds to leaving only one component of the field n fixed at the 
boundary, while by Dirichlet boundary condition it should be understood that the field 
n is fixed at the boundary, n(x,t)\ x=0 = n . 

Since the S-matrixof the GN model is, up to a CDD factor, the same as the 5-matrixof 
the O(N) nlcx model, the boundary Yang-Baxter equation (BYBE) for the two models 
is exactly the same, and so the reflection matrix of the GN model is given by the one of 
the nlcr model multiplied by the appropriate CDD factors^]. This indicates that there is 
a one-to-one correspondence between the integrable boundary conditions for the GN and 
nlcr model. Similarly to the bulk case, the different physics of these two models resides 
in these CDD prefactors. 

Let us summarize the results of the analysis of integrable boundary conditions for the 
GN model ||. The action of the boundary GN model is given by 

dx / dt C gn + dt C b , (3.1) 

-oo J— oo J— oo 

where Cb is the boundary action. As we have shown in the boundary lagrangian 

N i 

A = E^* > ( 3 - 2 ) 
i=i 1 

where ej = ±1, preserves the integrability of the GN model at the quantum level. The 
boundary condition derived from this action is 

1p+\x=0 = U=o • (3.3) 

Borrowing the terminology from the nla model, we will refer to bc's with all q = +1 or all 
€i = —1 respectively as Neumann and Dirichlet bc's. Therefore, up to index reshuffling, 
we have N + 1 inequivalent boundary conditions, which have to correspond to different 
solutions of the boundary Yang-Baxter equation. Due to the fact that the boundary 
interaction does not involve any flavor- changing terms, we should be able to find diagonal 
solutions for the BYBE, which will be done in section 4. In this section, following Ghoshal 
analysis, we exhibit solutions for the free and fixed boundary conditions, which serves as 
a warm-up for the more general case. 

2 This will give us a minimal solution, without any extra poles in the physical strip. 
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3.1 Free Boundary Condition 

As it follows from our discussion above, we should assume that the reflection amplitude 
for the elementary fermions is given by 

Ri(9)=R(9)6i. (3.4) 

Physically it means that we can not change the 'flavor' of the fermion by scattering it off 
the boundary and that the amplitude of scattering does not depend on the index i. For 
the nl<7 model this ansatz corresponds to no interactions on the boundary Q. In the case 
of GN model, it corresponds both to the bc's in equation with all equal to "+" or 
all €i equal to "— " , the difference between the two cases lying in a CDD factor. Due to 
the similarity between GN and nlcr model S'-matrices, R{9) can be written as 

R(8) = f(0)Ro(9), (3.5) 

where Rq(6) is the reflection amplitude for the nlcr model with free be computed by 
Ghoshal in [|: 

p (n\ _ V 2 ' 4tt 2tt / V ~ 2ir > V 4 _ 4?r ~ 2n > V 4 2tt / In a\ 

°^ ' r(i _j_ A _i_ M.\ v(l — r(3 _i_ A it\ rfl _i_ iL\ ' ^ ' 

1 \2 ' 4tt 2tt/ 1 V 1 2W 1 U 4tt 2tt/ 1 U 2tt/ 

The prefactor f(6) is fixed by unitarity and boundary crossing-unitarity, which generally 
read 

Rt(e)R{(-6) = 5l , K i \9) = Sp j ,(29)K i 'i'(-9), (3.7) 

where K l ^{9) = C u ' R\,{^ — 9), and C l i is the charge conjugation matrix. Equations |3.7| 
imply that f(9) should satisfy 

/(«)/(-«) = 1, (3.8) 

* = gj^§±i^ /g + ,). 

2 sinn 20 — i sm A 2 



The minimal solution of |3.8| , |3.9| can be found by elementary methods. In fact, there are 
two minimal solutions, with rather different physical properties. They are given by 

m ~ m sinh^-f) sinh^ + f + f) ' (3 - 10) 

and 



iX 



m = m sinh !^ + | sinh fr i 1 2 : , (3.11) 

1 ; sinh§(0-f) sinh±(# + f -f ^ 



differing by the CDD factor 



Fcdd(0) " tanh(f - | - f ) ' (3 ' 12) 
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The first solution exhibits only a simple pole in the physical strip, at 9 = i\/2, corre- 
sponding to the bound-state pole at 9 = iX in the bulk S-matrix. The second solution 
exhibits an additional simple pole at 9 = m /2 — i\/2, meaning that the boundary state 
for this solution contains a zero-rapidity single-particle contribution by the particle B j|] • 
In equations |3.10| , |3.11| $(0) is a prefactor satisfying 



?7T 77T 

*H0 = 1 , <$>{--9) = -<S>{- + 9). (3.13) 

The above equations are exactly those for the reflection matrix of the Ising model 
Since we do not expect any free parameters in our reflection matrices (the boundary 
term in the lagrangian |3.2| having no free parameters), we pick up the minimal solutions 
corresponding precisely to the boundary conditions ip + \ x=0 = ±ifj_\ x=0 for Ising fermions: 

, (m 9\ , (m 9\ 

$ + = - % coth f — — - J , $_=itanhl---j . (3.14) 



Therefore, we propose that R{9) in equation [3.5| with $ + corresponds to Neumann be 
and with $_ to Dirichlet be. 

As we noticed before, the reflection amplitude R(9) has a pole at 9 — y , corresponding 
to the bound-state pole at 9 = iX in the bulk S'-matrix. Since the particles B and B^ can 
be interpreted as bound-states of AiAj, we can use the boundary-bootstrap equation to 
compute the reflection amplitudes for them (see also fLl|), assuming that the boundary 



has no structure and consequently that the only non vanishing reflection factors are R^{9) 
and R B ^.{9). Recall that in general, if the particle can be interpreted as a bound-state 
of A ai A a2 (corresponding to a pole in the bulk scattering amplitude at 9 = iuj ), the 
boundary S'-matrixelements for Ab can be obtained by taking the appropriate residue at 
the bound-state pole of the two-particle boundary S-matrixi?"^^^ q^) M: 



where u = n — u and are the three-particle on-shell couplings defined by the residue 
of the bulk S'-matrixat the bound state pole: 

9 — iu b 

S C aZ(0)^ - fb aia2 fr 2 - (3-16) 

^ J 0,10,2 

Notice that the fused reflection amplitude is manifestly unitary, and we need only 
check boundary crossing-unitarity. A straightforward bootstrap computation gives 

R B 4(0) = --^-R(9„) R(9 + ))a 2 (29) , (3.18) 
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where 9 ± = 9±%. 

We can check the consistency of the bootstrap computation by verifying that the 
appropriate boundary crossing- unit arity equation is satisfied by these reflection ampli- 
tudes. It can be written easily, but one should be careful with factors coming from charge 
conjugation (see the appendix in The final result is 

K BB {6) = K BB (-9)S BB (29) + N{N ~ 1] K B ^(-9)S BB Bij (26) (3.19) 

where K BB (6) = flf (f - 9) and K BijBij (9) = -2i?gj(f - 9). By using the bulk 
amplitudes listed in the appendix, equation |3.19| can be easily shown to be satisfied. 



3.2 Fixed Boundary Condition 

Now let us consider the case where the first N — l fermions satisfy the "+" be and the N- 
th fermion satisfies the "— " be. In terms of reflection matrices, this situation is described 
by the ansatz 

Rt(9) ee R x {&), i = l,...,N-l, 

R%(6) ee R 2 (9), (3.20) 

which is the same as the one for fixed be in the nla model considered by Ghoshal, and 
therefore we follow his analysis closely. The amplitudes |3.20| can be written as 



R 1 {9) = f{9) R\{9), 

R 2 (9) = f(9)R° 2 (9), (3.21) 



where f(9) is given by |3.10| and R® and R 2 are the amplitudes for the nla model. From 



the BYBE Ghoshal found 

_ Ri{9) _ in -29 

x{6) = rW)~ ' ( } 



and by solving unitarity and crossing-unitarity 

Rl(0) - 



M2 T 4ir 2-n-J ^ ~ 2w > L \ 4 ~ 4vr ~ 2vr > L U 2tx > 



Rl(9) 



L \2 ~ 4tt ~ 2-k> ^ 2-k> \4 ~ 4vr 2tt > V4 ' 2v ) 

1 V2 4?r 2ti-/ 1 V 1 2tt> 1 U ±_ 4tt 2tt > 1 V 4 2tt / 

Wl j A I id_\ -p/i id_\ TV I j A id_\ -p/ 1 , id ■* 

1 V2 "T" 4tt "t" 2tt/ 1 V 1 2vr/ 1 U 4tt 2W 1 V 4 2tt / 



Using these reflection amplitudes we can compute, as before, the reflection factors for the 
two-particle bound-states. Notice that in this case we can obtain B by fusing fermions 
satisfying "+" be: 

R B (9) = i2x(0_) Rt(9 + ) [(N - 1) a x {29) + a 2 (29) + a 3 (29)} + R x (0-) R^a^) ; 

(3.23) 
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but we can also write B as the fusion of the iV-th fermion, which satisfies "— " be: 

R B (9) = R 2 (9_) R 2 {9 + ) [ax (29) + a 2 (29) + a 3 (29)] + R 2 (9-) Ri(9 + ) (N - 1) a x {29) . 

(3.24) 

These two expressions for R B (9) have to be equal, and this provides a non-trivial consis- 
tency condition for the boundary-bootstrap: 

X{6J) X(9 + ) [(N - 1) at(29) + a 2 (29) + a 3 (29)} + X(6-) ^(20) = 

= [ax{29) + a 2 (29) + a 3 {29)) + (N — 1) a t {29) X{9+) . (3.25) 



As we have checked, equation |3.25| turns out to be an identity. Notice that in this 
equation all information needed is the ratio R\/R 2 , which is fixed by the BYBE, and the 
ratio between bulk S'-matrixelements (CDD factors cancel out) which is fixed by the bulk 
Yang-Baxter equation; in other words, it depends only on the O(N) structure and not 
at all on CDD factors. This is quite surprising, since the consistency check is meaningful 
only if bound-states exist, which instead depends crucially on the presence of the CDD 
factor. 

The explicit expression for R^{9) is 

^- ^m^^ ^-^ 9 ^- (3 - 26) 

Similar bootstrap computations give the reflection amplitudes for the antisymmetric ten- 
sor components, 

iffi(g) = - {tX ^ e 26) Ri(0-)Ri(e + )a 2 (29) , i,j*N, (3.27) 

{i\ + 29){m - i\ + 29) 
R B f N {e) = m in-i\-20) fr( g -W g +M2g) • (3-28) 

These amplitudes satisfy unitarity automatically. In this case boundary crossing-unitarity 
reads 

K BB {9) = K BB {-9)S B B B B {29) + (A r ~ 1) 2 (A T ~ 2] - K B ^(-9)S BB Bij (29) + 

+ (N- l)K B ^(-9)S BB Bm (29) . (3.29) 

and we have checked that it is indeed satisfied. 



4 General Boundary Conditions 



For the general boundary condition |3.3| let us assume that for i = 1, . . . , M, 6j = +1 
and for i — M + 1, . . . , N, 6j = —1, which is, up to index reshuffling, the most general 
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boundary condition we have to consider. Since, as before, the boundary action has no 
term involving flavor change, we have to assume that the reflection matrix is diagonal. 
Therefore we have the following ansatz for the non-vanishing elements of the reflection 
matrix^ 

Rl(6) = R 1 (e), R a a (8) = R 2 (8) . (4.1) 

Similarly as in previous case, the BYBE fixes the ratio Ri(9)/R 2 (9). Consider first the 



factorization of the two-particle reflection process | Ai(#) Ai(9'))i 
it gives the following equation: 



A(-e') A a (-8)) out ; 

(4.2) 



a 2 (M Ri(8i) a 3 {9 12 ) R 1 {9 2 ) + <r 3 {0 12 ) R 2 (6 1 ) a 2 {9 12 ) R 1 {9 2 ) = 

= R2{6 2 ) R 1 {9 1 ) a 2 (9 12 ) R^Bx) a 3 (9 12 ) + R 2 (6 2 ) a 3 (9 12 ) R 2 {9 X ) a 2 (6 12 ) , 

where 9% 2 = Q\ + 9 2 . By dividing this expression by R 2 {9i) R 2 {9 2 ) o~ 2 {9i 2 ) a 2 {9i 2 ) and 
taking the limit 9\ — > 9 2 we get a differential equation for X{9) = Ri(9)/R 2 (9), 



d X(9) - X2{9) ~ 1 
d9 X{9) ~ 29 



whose solutions are 



X{9) 



C-9 
C + 8 



and 



X{9) 



(4.3) 



(4.4) 



C being an arbitrary integration constant. The solution X(9) = 1 corresponds to Neu- 
mann or Dirichlet boundary conditions, since in these cases R\{9) = R 2 (9), which we 
have analyzed in the previous sectionFj. 



Consider now the factorization of the process \A a {9) A a {9')) ir , 
It gives the following equation: 



\A a {-9')A a {-9)) 



out- 



[Ma 1 {9 12 )+a 2 {9 12 ) + a 3 {9 12 )}R 1 {9 1 )a 1 {9 12 R 2 {9 2 ) + 

+ <xi(M R 2 (0i) \(N - M)ax(9 12 ) + a 2 (8 12 ) + a 3 (9 12 ) 



R 2 



Ri(9 2 ) Max 



7 12 



0-2W12 



0-3(^12 



Rx{9) *i(9 12 ) + 



+ R 1 (9 2 ) <J\{9\2) R2{9 1 ) [(N - M)a 1 {9 12 ) + a 2 (9 12 ) + a 3 (9 12 )) . (4.5) 



If we plug X{9) in this expression, the final, compact result is that C is fixed to be 

m N - 2M 



C 



(4.6) 



2 N -2 

By analyzing the BYBE for the other processes we find that [4.6| is the unique consistent 



solution. Notice that by taking M = N — 1 in equation 4.6 , we obtain C = in/2 and, for 

3 From now on we will use letters in the middle of the alphabet (i, j, ...) for fcrmions labeled from 1 
to M, and letters in the beginning of the alphabet (a, b, ...) for fermions labeled from M + 1 to N. 
4 Recall that the Neumann and Dirichlet cases will only differ by overall CDD factors. 



8 



the nlcr , we recover Ghoshal's results for fixed be. All that is left to do is to compute 
the prefactors for the reflection amplitudes, by using unitarity and boundary crossing- 
unitarity |3.7| , which will fix R\(9) and R 2 {9) up to CDD factors. These conditions read 
explicitly 



R 1 {e)Ri{-t 



1 

iX 



29 



i\ + 29 



R 2 (9)R 2 (-9) = 1 
iX(N-M- 1) - 29 
iX(N - M - 1) + 29 



<rj(20) #1 



ATT 



+ 



(4.7) 
(4.8) 



where cri(9) = Na\(9) + a 2 {9) + 03 (6>). The minimal solution of equations Hl| is given by 



Tfi + — + — 
R 1 (9) = -f(9)R (9) U !; - 



&-&)ni + &(N-M-i) 



+ — 



X 



M-l)-f) 



r(| + ^(iv-M-i)- 



i§_- 
2w 



2tt' 



(4.9) 



and 



rfi + — + — 

R 2 (9) = -f(9)Ra(9) U - ^ 



U T 4i 27r^ U 



A(iV-M-l) + f 



x 



1 V 4 ~ 4vr 27T / V 4 ~ 47T ~ 27T > 

n-i+&N-M 



A.(N-M-l)-£ 



(4.10) 

The same amplitudes, but 



+ 

where /(#) is given in equation |3.10| and _R i n equation 
with / Rq replaced by R , apply to the nlcx model. 

Notice that R 2 has a pole at 9 = — ^(iV — 2M), which is in the physical strip for 
N/2 < M < N — 1. As argued by Ghoshal, this pole signals the presence of one-particle 
contributions in the boundary state. Since upon the substitution M — > iV — M, the ratio 
X{9) -> l/X{9) and i2j and R 2 get interchanged, for 1 < M < N/2 it will be R x to 
exhibit this pole in the physical strip. 

It is interesting to notice also that for M = N, R\ reduces to the second minimal 



solution for the free boundary condition, equation 3.1'J 



The amplitudes for the two-particle bound-states can be computed by boundary- 
bootstrap, and they read explicitly 

[iX + 29){2C -iX + 29) 



R%{9) 



(29)(2C 
[iX + 29) 



-R 1 (e.)R 1 (e + )a 2 (29) , 



29 

{iX + 29)(2C 



iX - 29) 
R 1 (9_)R 1 (9 + )a 2 (29), 

tX + 29)(2C + iX + 29) 



(29)(2C -iX- 29){2C + iX- 29) 
{iX + 29)(2C - iX + 29) 



R X (Q-)R 



H (7 + 



(29)(2C - iX - 29) 



-R 1 (9^)R 1 (9 + )a 2 (29) 



(4.11) 
(4.12) 
)<r 2 (20),(4.13) 
(4.14) 



In this more general case, the appropriate bootstrap consistency condition corresponding 
to the generalization of equation |3.25| is 



X(0_)A:(0 + )[M<7i(20) + a 2 (20) + o 3 {26)] + X{6„){N - M)a 1 {26) = 

= [(N - M)a x {29) + a 2 {29) + a 3 (20)] + Ma 1 (29)X(9 + ) . (4.15) 



It is very easy to see that if M = N — 1 this reduces to equation |3.25|. As before we should 



stress that bootstrap consistency requires only information obtained from the BYBE, via 
the ratio X(9), and the fact that the elementary fermions form isoscalar bound-states. 
Finally, we have also explicitly checked that the reflection amplitudes listed above satisfy 
the appropriate boundary crossing-unitarity. 



5 Conclusions 

In this paper we have computed the minimal boundary S-matrixfor the Gross-Neveu 
model, extending Ghoshal's analysis of the nonlinear a model. We found general (diag- 
onal) solutions for the boundary Yang-Baxter equation for both models and connected 
them to the boundary conditions proposed recently in our paper || . We also proved that 
the solutions presented in this paper are consistent with the boundary-bootstrap. This 
seems to indicate that the boundary contains a lot of information of the bulk theory. It 
would be interesting to investigate how much we can learn about bulk integrable models 
starting from the boundary. As a natural follow-up to this work it would be interesting 
to study the boundary Yang-Baxter equation in general and see if there are non-diagonal 
solutions, that is flavor- changing scattering off the boundary. There should be possible, 
then to find associated microscopic boundary conditions for the elementary fields. 
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A Appendix 



In this appendix we list the S'-matrixelements for scattering processes involving two- 
particle bound-states (the isoscalar and the first antisymmetric tensor), that we need to 
write down the boundary crossing-unitarity equation |3.19| , |3.29| . These amplitudes are 
obtained from those of elementary fermions by fusion. 

Recall that if the scattering amplitude of two particles A a and Ab has a simple pole 
in the physical strip at, say, 9 = iu c ab) with residue given by : 



Si v {9) 



c fa'b' 



1 fab ft 



9 — iu c „ 



(A.l) 



then the scattering of the bound-state particle A c with all other particles in the theory 
can be obtained by the bootstrap equation |L(| 



fa b S C cf(9) = f c aibl S b h t{9 - iul) S a X(9 



~ia\d' , 



(A.2) 



where u = it — u. In the case of GN model, we get 

9(i7t - 9) [A(tt - 3A) + 9(iir - 9)] - 2ttA 2 (A - tt)) 



S 



BB i 
BB 1 

qBij Bij 

^BB 



s 



BB 



QBBij 
^BB tJ 1 



9{m - 9){9 - i\)(iir - 9 - iX) 
2i(N-A)X 3 



a 2 (9 + )a 2 (L)a 2 2 (9) 



(in - 9)(9 - i\)(m - 9 - iX) 

i{N-4)X 3 

2(m - 9)(9 - iX){m - 9 - iX) 
6{m -9) + A(tt - 3A) 



<7 2 (0 + )<7 2 (0_)<7 2 2 (0) 



a 2 (9 + )a 2 (L)a 2 2 (9) 



- 9 - iX) 
-i{N-A)X 3 



9(9-iX)(in 



iX) 



a 2 (9 + )a 2 (9_)ai(9) : 
a2(9 + )°2(9-)°l(Q) 



where 9± = 9 ± iX. Unitarity follows directly from the fusion procedure and crossing 
symmetry is satisfied with charge conjugation matrix elements C BB = 1 and C BijBij = 
-2. 



11 



References 



[1] S. Ghoshal and A. Zamolodchikov, Boundary S-matrixand Boundary State in Two- 
Dimensional Integrable Quantum Field Theory, Int. Jour. Mod. Phys. A 9, (1994) 
3841-3886; |hep-th/ 9306001 . 

[2] A. Fring and R. Koberle, Factorized Scattering in the Presence of Reflecting Bound- 
aries, Nucl. Phys. B 421, (1994) 159-172; [Eep-th/ 9304141 . 



[3] D. Gross and A. Neveu, Dynamical Symmetry Breaking in Asymptotically Free Field 
Theories, Phys. Rev. D10, (1974) 3235-3253. 

[4] S. Ghoshal, Boundary S-matrixof the O(N)- Symmetric Non-Linear Sigma Model, 
Phys. Lett. B334, (1994) 363-368; |hep-th/9401008| . 

[5] M. Moriconi and A. De Martino, Quantum Integrability of Certain Boundary Con- 
ditions; |hep-th/9809T78 . 



[6] R.F. Dashen, B. Hasslacher and A. Neveu, Semiclassical bound-states in an asymp- 
totically free theory, Phys. Rev. D12, (1975) 2443-2458. 

[7] A. B. Zamolodchikov and Al. B. Zamolodchikov, Exact S Matrix of Gross-Neveu 
"Elementary" Fermions, Phys. Lett. B72, (1978) 481-483. 

M. Karowski and H.J. Thun, Complete S-matrixof the 0(2 N) Gross-Neveu Model, 
Nucl. Phys. B190, (1981) 61-92. 

[9] E. Witten, Some Properties of the (ipip) 2 Model in Two Dimensions, Nucl. Phys. 
B142, (1978) 285-300. 

[10] A. Zamolodchikov and Al. Zamolodchikov, Factorized S-Matrices in Two Dimen- 
sions as the Exact Solutions of Certain Relativistic Quantum Field Theory Models, 
Ann. Phys. 120,(1979) 253-291 . 

[11] N. MacKay, Fusion of SO(N) Reflection Matrices, Phys. Lett. B357, (1995) 89-93; 
|hep-th/9505124 . 



[12] A. De Martino and M. Moriconi, Tricritical Ising Model with a Boundary, Nucl. 
Phys. B528, (1998) 577-594; |hep-th/9803136| . 



12 



